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The pp — > pmr+ reaction is a channel with the largest total cross section for pp collision 
in COSY/CSR energy region. In this work, we investigate individual contributions from 
various N* and A* resonances with mass up to about 2 GeV for the pp — ► pmr + reaction. 
We extend a resonance model, which can reproduce the observed total cross section quite 
well, to give theoretical predictions of various differential cross sections for the present 
reaction at T p = 2.88 GeV. It could serve as a reference for identifying new physics in 
the future experiments at HIRFL-CSR. 



1. Introduction 

The study of excited nucleon states is very important for understanding the inter- 
nal structure of nucleon and the strong interaction in the nonperturbative QCD 
domain In the early years, our investigation on the N* and A* baryon spec- 
troscop y w as mainly based on ttN experiments, which made observations unsatis- 
factory^^. An outstanding problem is that, in many of its forms, the quark model 
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predicts a large amount of "m issing " N* and A* states around 2 GeV/c 2 , which 
have not to date been observed lalslal Therefore, it is of necessity to search for these 
"missing" N* and A* states from other production processes. Moreover, even for 
those well-established resonance states, properties like mass, width and branching 
ratios still suffer large experimental uncertainties 1 ^, which also need further studies 
in more other production processes. Here we propose to look for "missing" N* and 
A* resonances in pp — > pnn + reaction. At COSY/Julich, Experiments for studying 
N* and A* resonances through pp collisions are being carried out, but there is a 
lack of a good An detector for complete measurement of diverse differential cross 
sections. At present, a heavy ion cooler-storage ring HIRFL-CSR — an accelerator 
system of the same beam energy region with maximum incoming-proton kinetic en- 
ergy up to 2.88 GeV^, has already been completed at Lanzhou. With its scheduled 
4-7T hadron detector it will have a special advantage for studying excited nucleon 
states through pp collisions. 

Recently, BES collaboration has produced quite a few novel findings on N* reso- 
nances by using various N* production processes from J/ij) or ip' decays^ ^^^^ . 
Measurements of the J/ip — > pn^n + c.c. decay by BES collaboration showed two 
new, clear N* peaks in the pir invariant mass spectrum around 1360 MeV/c 2 and 
2065 MeV/c 2 , respectively 121 Of them the former one was identified as the first 
direct observation of the A* (1440) peak in the ttN invariant mass spectrum, which 
was confirmed by the CELSIUS- WASA Collaboration in their nir + invariant mass 
spectrum of pp — > pnn + reaction For the latter one, it is very likely to be 
a long-sought missing N* peak aroun d 2 GeV/c 2 . However, similar searches for 
it in ip' decays are inconclusive EHH . Therefore, it is of necessity to look for the 
new N* resonance i n ot her reaction processes, such as the pp — > pnir + reaction. 
Furthermore, in Ref. 031, the authors found that the A*++(1620) resonance gives 
an overwhelmingly large contribution in the pp — > nK + Tj + reaction by t-channel 
p + exchange. If so, it is also expected to make a significant contrib utio n in the 
pp — > pnn + reaction, as can be checked in the present work. In Ref. ^| we have 
studied the pp — > pnir + reaction for beam energies below 1.3 GeV. Here we extend 
the study of this reaction to higher energies and investigate individual contributions 
from various N* and A* resonances with mass up to 2 GeV/c 2 for this reaction. We 
extend a resonance model, which can describe the experimental data of the total 
cross section for beam energies ranging from 0.8 GeV to 3.0 GeV quite well, to 
give theoretical prediction of various differential cross sections for the pp — > pnir + 
reaction at T p = 2.88 GeV. It can be used for the subsequent comparison with 
the experimental results at COSY and HIRFL-CSR. Meanwhile, it could serve as a 
reference for the construction of the scheduled 4ir hadron detector at HIRFL-CSR. 

2. Formalism and ingredients 

We study the pp — > pnir + reaction within an effective Lagrangian approach. In 
our model, all the mesons, baryons and resonances are treated as fundamental fields. 
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All the basic Feynman diagrams involved in our calculation for this reaction are 
depicted in Fig. [T] In view of overall system invariant mass about 3 GeV for T p = 
2.88 GeV, we have checked contributions from all the well-established N* and 
A* resonances (overall status 3 or 4 stars) below 2 GeV/c 2 , but only present the 
results of the relatively significant ones in next section. Meanwhile, we investigate 
the contribution from iV*(2065) resonance for the present reaction. Explicitly, we 
list in Table [T]all the N*, A* resonances and the meson exchanges considered in 
our present calculation. 

For TV* (2065), according to results in Ref.^, its spin-parity is limited to be l/2 + 
and 3/2 + , and it is more likely that both are needed. In the quark model there are 
predictions for the existence of N * resonances with spin-parity l/2 + and 3/2 + 
between 2.0 and 2.1 GeV/c 2 [SEE! Since the spin-parity of the new resonance(s) 
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was not well determined, we assume that this peak consists of exactly those reso- 
nances with JP = l/2 + ,3/2 + predicted in Ref. 5 , which are A*(1975)(J P = l/2 + ), 
A* (2030) (J p = 3/2+) and A*(2065)(J P = l/2 + ). Among them A* (2065) has 
much stronger coupling to irN than the other two, which is in accord with BES 
results. As in Ref.H^I here we also treat the observed A''* (2065) peak as an effective 
A*(2065)(J P = l/2 + ) resonance which represents all contributions of the three res- 
onances AT* (1975) (J* = 1/2+), A*(2030)(J P = 3/2+) and A*(2065)(J P = 1/2+). 
In so doing, the coupling constant g 2 T T NN*{2065)/^ 7r is scaled by a factor of 1.122. 
See Ref. for details of this effective treatment. Of course, we have used the BES 
observed values for the mass and width of A* (2065) to determine its relevant cou- 
pling constant, see Table [TJ In view of scanty information for its decay branching 
fractions, here we regard A-7T as the dominant decay mode of A* (2065), whereas 
the A* (2065) peak in invariant mass M p7T - spectrum is so strong and highly sig- 
nificant. So, in our calculations we have taken an artificial branching ratio up to 1 
for Att decay mode. 

The effective Lagrangian densities involved for describing the meson- AA ver- 
tices are: 

£~xNN = ~ ^ NN UNlhlfif ■ dPtpnUN, (1) 



C V NN = ~igr ] NNUNl5i l i 1 UN, (2) 
£<jNN = g<jNNUNlpaUN, (3) 



CpNN = -9 p nnun(Hp + Op V d v )T ■ ip^UN- (4) 

At each vertex a relevant off-shell form factor is used. In our computati on, we 
take the same form factors as used in the well-known Bonn potential model ^* 

= ( A Jr ? 2 V (5) 

with n=l for 7r, 77 and a mesons and n=2 for p meson, ku, m-M and Am are 
the 4-momenta, mass and cut-off parameter for the exchanged meson (M), respec- 
tively. T he coupling con stants and the cut-off parameters are taken as the follow- 
ing ones U4 | 15 | 17 | 18 M 9 l NN /4n = 14.4, g$ NN /4* = 0.4, A, = A„ = 1.3 GeV, 
gl NN /kn = 5.69, A CT = 2.0 GeV, g 2 pNN /A-K = 0.9, A p = 1.85 GeV, and k = 6.1. Note 
that the constant g^NN is related to / w nn of Eq.fl} as g-^NN = {firNN /m w )2m,N&^. 

To calculate the amplitudes of diagrams in Fig. [1] within the resonance model , we 
also need to know interaction vertices involving A* and A* resonances. In Ref. ^H, 
a Lorentz covariant orbital-spin scheme for A* AM couplings has been described 
in detail, which can be easily extended to describe all the couplings that appear 
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in the Feynman diagrams depicted in Fig. Q] By using that scheme, we can easily 
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obtain the effective couplings: 



£-ttNA 


(1232) 




.9A(1232)7V7r«A r <9 A V?r'UA(1232) M + h -C, 


(6) 


£-rtNN* 


(1440) 


= 


9n* 


(1440)Af7r"Af757M^ ? /'7rWAf«(1440) + h.C., 


(7) 


£<jNN* 


(1440) 
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9n* 


(1440)AT<t%^(7%« (1440) +h.C, 


(8) 




(1520) 
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9n* 


(1520)NirUN757tJ.P%Pn' l PwUN* (1520)i/ + h.C., 


(9) 


r 

'"pNN* 


(1520) 




CI ivr* 

yN* 


1 -\ rnn\ jvt 11 TVrl/?^ 11 AT* I -\ r on\ — I— ll C 

(152Q)Np Uj N y p a N*{152Q)fi i 11.^., 


(10) 


r 

>~"kNN* 


(1535) 




CI ivr* 

yN* 


(1535) Ntt u, iV tyir * (1535) ~ n.^.j 


fill 


r 

'-"qNN* 


(1535) 




yN* 
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(12) 


£pNN* 


(1535) 


= 


9n* 


(1535)NpUN-f5(lp ^2 )1pp U N * (1535) + h.C., 


(13) 


£-rtNA* 


(1600) 


= 


9a* 


1600)iV7rWlV<9' i 'i/VUA*(1600)M + h.C., 


(14) 


£-ttNA* 


(1620) 
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9a* 


1620) Nt^N^-rU^ (1620) + h.C, 


(15) 
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(16) 
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(17) 


n 
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(18) 
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(1650) 
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(19) 


C-kNN' 


(1675) 


= 


9n* 


(1675)Af7rWATP^K 1 /'7 r ' u Ar*(1675)^ + h-C, 


(20) 


£-irNN* 


(1680) 
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9n* 


(1680)JV7rWAr757^KP7rV , 7rMjV*(1680)i.A + h.C, 


(21) 
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For the relevant vertices involving N* and A* resonances, the off-shell form 
factors are adopted as follows: 

where n=l for all the resonances except for n=2 for A(1232). All the coupling 
constants and cut-off parameters used in the present paper are listed in Table [TJ 
In addition, we also introduce form factors for the off-shell baryon resonances as in 

F «® = A* + {*-m> R r (38) 

with A= 0.8 GeV. 

The propagators can be written as 

4+m R 
q 2 - m\ + im R T R 



G{9) = S Y . (39) 

— m 1 ■ ■ ■ 1 



for the spin-i resonances, 



gm = - — f M r (g) r nd) 

q z -m R + im R T R 



with 



11 2 

P»v(q) = -{4 + m R )[g^ - -7^7^ - ^ {l^u - 1^) ~ ^— r^*,], 



6m R 6m 



for the spin-| resonances, and 



with 



(41) 



<W(<?) = rr r ( 42 ) 

g z - m R + im R l R 



P^aff{q) = - {4 + m R )\^{g^ a g u p + g^pg ua ) ~^g^ v g a p (43) 
+ (7^7^/3 + lvlfi~g m + %1l39va + >7a3ju/3)L ( 44 ) 

- , s . q»qv - . 

SV (g) = ~9tiy + — 2~, 7m = ~7m h — r- I 45 ) 
m R m R 

for the spin-| resonances. 

After the effective Lagrangians, coupling constants and propagators fixed, the 
amplitudes for various diagrams can be written down straightforwardly by following 
the Feynman rules. And the total amplitude is just their simple sum. Here we give 
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Table 1. Relevant parameters of N* and A* included in our calculations. The widths and branching 
ratios are taken from PDG 6 and the cut-off parameters are from Refs. 14 ' 15,18 ' 25 . Here the g 2 /4tt 
for A r *(2065)A f 7r vertex has already been scaled by a factor of 1.122. 

Resonance Width/GcV Decay mode Branching ratio g 2 /^ Cut-off/GeV 
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0.012 
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0.042 
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0.8 
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Ntt 
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0.18 
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1.3 



explicitly the individual amplitude corresponding to 7V*(1440)(7r° exchange), as an 
example, 

X(iV*(1440),7r ) = V2^^g 2 N , Nn u n (p n ,s n ) 75 ^G N * ^ {q) lh }i w u 1 (p 1 ,s 1 ) 



in — " k'i ill 

U3(P3, S3)75 /^7rW2(p2, s 2 ) + (exchange term with pi <-> p 2 ), (4( 
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where u n (p n , s n ), ^3(^3, S3), u±(px, s±), 1*2(^2, S2) denote the spin wave functions of 
the outgoing neutron, proton in the final state and two initial protons, respectively. 
p n and fcjr are the 4-momenta of the outgoing and the exchanged pion mesons, p\ 
and P2 represent the 4-momenta of the two initial protons. The coupling constant 
appearing herein can be determined from the experimentally observed partial decay 
width of iV*(1440) resonance as follows, 

^9 2 n-(uao)NttP C n 1 ,m^{E N - m N ) 



with 



IV* (1440)^^ 



-in 



m N*(1440) 



2{p' 



Pn 



\ 



( m W«(1440) _ ( m N + TO ^) 2 )( m W-(1440) ~ ( m N ~ m ^) 2 ) 



(48) 



En — 



(49) 



All the other coupling constants can be obtained similarly. See Refs. | 14 | 15 | f or 
details. 

Then the calculation of the cross section a{pp — > pnir + ) is straightforward, 
, , + , 1 ™>l 1 A ^ "V^ d 3 pn m n d 3 p n 4 

da(pp^ pnn^) = -— \M\ — — d {pi+P2~P3-Pn-p n ) 



Si,Sf 



E 3 2E V E n 



with the flux factor 



F=(27r) 5 v /(p 1 -p 2 ) 2 -m4. 



(50) 



(51) 



The factors 1/4 and Yl s . Sf emerge for the simple reason that the polarization of 
initial and final particles is not considered. 



3. Numerical results and discussion 

With the formalism and ingredients discussed in the former section, we com- 
puted the total cross section versus the kinetic energy of the proton beam (Tp) 
for the pp — > pnir + reaction by using a Monte Carlo multi-particle phase space 
integration program. The results for Tp ranging from 0.8 to 3.0 GeV are shown in 
Fig. [5] along with experimental data for comparison. 

As one can see from Fig.[2j the experimental data of total cross section are repro- 
duced reasonably well by our theoretical calculations over the entire energy range. 
Note that we have considered the interference terms between the direct amplitudes 
(diagram a,c in Fig. [1]) and the corresponding exchange amplitudes (diagram b,d 
in Fig. [1]) in our calculations. However,the interference terms between different 
resonance-excitation processes and between various meson-exchange diagrams are 
ignored. We also show contributions of various components which are large and not 
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1.8 2.0 2.2 2.4 

T„(GeV) 



2 6 2.8 3.0 



Fig. 2. Total cross section and contributions from various TV* (left) and A* (right) resonances 
as a function of Tp for the pp — > pnn + reaction with the solid line as the simple incoherent 
sum of all contributions, compared with data 26 . Left: the dashed, dotted, dot-dashed, dot-dot- 
dashed, short-dashed and short-dotted lines represent individual contributions from TV* (1440), 
TV* (1520), TV* (1650), N* (1675), TV* (1680), and TV* (2065), respectively. Right: the dashed, dot- 
ted, dot-dashed, dot-dot-dashed, and short-dashed lines represent individual contributions from 
A(1232), A* (1600), A* (1620), A* (1700), and A* (1905), respectively. 



negligible there, N* contributions in Fig. [5] (left) and A* contributions in Fig. [2] 
(right), respectively. Individual contributions from iV*(1440), iV*(1520), V*(1650), 
A*(1675), iV*(1680), an d iV*(2065) are presented in Fig. [3 (left) by dashed, dot- 
ted, dot-dashed, dot-dot-dashed, short-dashed and short-dotted lines, respectively. 
And contributions from A(1232), A*(1600), A*(1620), A*(1700), and A*(1905) 
are shown in Fig. [2] (right) by dashed, dotted, dot-dashed, dot-dot-dashed, and 
short-dashed lines, respectively. One can find that contributions from A(1232) and 
N* (1440) are still dominant in present energy region and the contribution from 
iV*(1680) becomes significant for kinetic energy above 2.0 GeV. We also give our 
predictions of invariant mass spectra and Dalitz plot in Fig. [3] and the momentum 
and angular distributions of the final charged particles in Fig. [4] for pp — > pnir + 
reaction at T P = 2.88 GeV. 

The pp — > pmr + reaction is a channel with the largest total cross section for 
pp collision in the present energy region. Since the kinetic energy Tp of the proton 



beam at HIRFL-CSR can reach 2.88 GeV with luminosity above 10 cm s 



10 



the event rate will be so large that it is easy to collect enough events and to get 
high statistics data for this channel at HIRFL-CSR. The theoretical prediction of 
various observables given in this work would provide us more knowledge on the 
relevant physics when compared to the future experimental results. Meanwhile, a 
good measurement of invariant mass spectra and Dalitz plot will play the key role 
for identifying new resonances and determining their parameters. In this regard, 
the scheduled Air hadron detector at HIRFL-CSR will be particularly competent. 
As mentioned above, the spin-parity of the iV*(2065) peak was not well deter- 
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1.5 2 1 1.5 

M pn +(GeV) M nn +(GeV) 




M pn (GeV) M'„ + (GeV 2 ) 



Fig. 3. Theoretical prediction of invariant mass spectra and Dalitz plot for pp — > pnn^ reaction at 
Tp = 2.88 GeV. The dashed curves stand for pure phase space distribution while the solid curves 
include the interaction amplitudes. 



mined by the BES collaboration, therefore we have u sed an effective treatment for 
the TV* (2065) resonance(s) as the authors did in Ref. 1^1 This effective description 
would indeed be a very good approximation for the total cross section, but generally 
speaking, it might reproduce the differential cross sections not so well. However, 
due to the dominance of iV*(2065) (its magnitude is much stronger than the other 
two) among the three resonances predicted in Ref. hence even for the description 
of various differential cross sections, it would be acceptable. Of course, this issue 
still waits for an exact answer from future experimental results at HIRFL-CSR. 

To sum up, in this paper we investigate individual contributions from diverse 
N* and A* resonances up to 2 GeV/c 2 for the pp — > pnir + reaction. We extend 
a resonance model, which can describe the observed total cross section for beam 
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Fig. 4. The momentum and angular distributions of the final proton and charged pion for pp — > 
pnw + reaction at Tp = 2.88 GeV, compared with pure phase space distributions (dashed curves). 



energies ranging from 0.8 GeV to 3.0 GeV quite well, to give theoretical prediction 
of various differential cross sections for this reaction at T p = 2.88 GeV. It can be 
used for identifying new physics in the future experiments at HIRFL-CSR. It could 
also serve as a reference for the construction of the scheduled Air hadron detector 
at HIRFL-CSR, which is quite possible to offer more physical information and to 
help us understanding the relevant physics better. 
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